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ABSTRACT 


A  ten-level  primitive  equation  ocean  circulation  model 
is  used  to  investigate  the  formation  and  evolution  of 
large-scale  thermal  anomalies  observed  in  the  central  North 
Pacific  Ocean  during  the  fall  and  winter  of  1976.   Several 
initial  value  model  integrations  of  4  months  duration  are 
carried  out  in  order  to  help  explain  the  observed  anomaly 
development.   Initial  and  verifying  ocean  data  down  to  400  m 
depth  are  obtained  from  the  NORPAX  ships  of  opportunity  pro- 
gram.  Anomalous  atmospheric  wind  forcing  is  obtained  from 
Namias '  monthly  mean  sea-level  pressure  anomalies,  while 
climatological  heat  fluxes  are  used. 

The  skill  with  which  the  model  simulates  the  observed 
anomaly  evolution  in  the  different  experiments  is  estimated 
synoptically  and  measured  statistically  by  calculating  root 
mean  square  (RMS)  temperature  errors  and  SI  skill  scores. 
Analysis  indicates  that  anomalous  atmospheric  wind  forcing 
improves  the  model  predictions  in  the  upper  levels.   For 
this  particular  winter  case  using  climatological  heating, 
however,  knowledge  of  the  initial  anomalous  temperature  con- 
ditions does  not  improve  the  model  results.   The  model  skill 
at  upper  levels  exceeds  both  persistence  and  climatology 
(forecast  of  zero  anomaly)  while  at  the  lower  levels  it  is 
comparable  to  persistence  but  not  as  good  as  climatology. 


TABLE  OF  CONTENTS 

I.  INTRODUCTION  ------------------9 

II.  DATA  PROCESSING   ---------------  -12 

III.  MODELING  APPROACH   --------------  -17 

IV.  RESULTS  AND  CONCLUSIONS   -----------  -22 

LIST  OF  REFERENCES  -----------------  -45 

INITIAL  DISTRIBUTION  LIST   -------------  -47 


LIST  OF  TABLES 


Table  I.    Definition  of  terms  used  in  equations 

(1)  and  (2)  --------------  -41 

Table  II.   Ratios  of  root  mean  square  error  (RMSE) 
for  persistence  (PERS)  and  ocean  models 
(1,2,3,4  and  5)  to  RMSE  for  climatology 
for  October  1976,  November  1976  and 
January  1977  at  10,32,62,102,162,262, 
462,  900  and  1700  meters;  followed  by 
similarly  nondimensionalized  ratio  of 
SI  scores  for  same  times  and  levels  -  -  -  42 


6 


LIST  OF  FIGURES 

Fig.  1.   Geographical  area  occupied  by  the  idealized 

rectangular  ocean  model;  the  TRANSPAC  (ocean) 

data  area  and  the  method  of  extrapolating 

the  oceanic  temperature  data  outward  into 

the  data  void  regions  of  the  model  along  the 

line  AO --  31 

Fig.  2.   Anomalous  pressure  fields  in  the  central 

Pacific  for  September  1976  to  January  1977 

at  2  mb  intervals;  shaded  areas  represent 

negative  anomalous  pressure;  area  shown  is 

the  area  east  of  the  dateline  in  which  the 

temperature  data  are  given  in  figures  3 

through  8-------------------  32 

Fig.  3.   Observed  temperature  anomaly  fields  for 

September  1976,  November  1976  and  January 
1977  at  10,  102,  162  and  262  meters;  con- 
tour interval  at  10  meters  is  0.5°C;  con- 
tour interval  at  102,  162  and  262  meters 
is  0  .  2  ° C ;  shaded  areas  represent  negative 
anomalous  temperatures   ------------  33 

Fig.  4.  Same  as  figure  3  except  for  Model  1  results  -  -  34 

Fig.  5.  Same  as  figure  3  except  for  Model  2  results  -  -  35 

Fig.  6.  Same  as  figure  3  except  for  Model  3  results  -  -  36 

Fig.  7.  Same  as  figure  3  except  for  Model  4  results  -  -  37 

Fig.  8.  Same  as  figure  3  except  for  Model  5  results  -  -  38 

Fig.  9.   Plots  of  root  mean  square  error  (RMSE)  for 
persistence  (P)  and  ocean  models  (1,  2,  3, 
4  and  5)   nondimensionalized  by  the  RMSE 
for  climatology  as  a  function  of  time  in 
days  for  10,  102,  162  and  262  meters   -----  39 


ACKNOWLEDGMENTS 

The  author  wishes  to  express  his  appreciation  to 
Dr.  Robert  L.  Haney  for  his  guidance,  kind  understanding 
and  counsel  in  the  preparation  of  this  thesis  and  to 
Dr.  R.  J.  Renard  for  reading  the  thesis  and  making  many- 
valuable  comments. 

Special  thanks  are  extended  to  Drs.  R.  L.  Bernstein 
and  W.  B.  White  of  Scripps  Institution  for  their  graciously 
allowing  the  use  of  the  TRANSPAC  XBT  data.   Appreciation 
also  goes  to  Dr.  J.  Namias  of  Scripps  for  the  use  of  his 
sea  level  pressure  anomaly  data  for  September  1976  to 
January  1977. 


8 


I .    INTRODUCTION 

Early  investigations  by  Namias  and  others  indicated  that 
the  winds  are  an  important  source  of  kinetic  energy  for  the 
oceanic  circulation.   But  just  as  the  mean  atmospheric  cir- 
culation helps  drive  the  mean  large  scale  ocean  circulation, 
anomalous  atmospheric  forcing  can  be  assumed  to  drive  anoma- 
lous currents  in  the  ocean  (Malkus,  1962) .   This  premise  was, 
in  part,  the  motivating  force  for  the  present  model  study. 
The  effects  of  anomalous  atmospheric  forcing  on  the  oceans 
have  been  a  matter  of  intuitive  guessing  for  years,  but  few 
studies  have  attempted  to  assess  in  more  than  a  descriptive 
or  qualitative  manner  the  nature  or  magnitude  of  these 
effects.   A  notable  exception  is  the  recent  work  of  Davis 
(1976,  1977) . 

One  of  the  most  important  manifestations  of  the  effects 
of  anomalous  atmospheric  forcing  on  the  oceans  is  the  devia- 
tion of  the  ocean  thermal  structure  from  the  long  term  cli- 
matological  norm.   The  space  and  time  scales  of  these  devia- 
tions or  anomalies  are  known  in  a  very  general  sense,  but 
the  exact  generating  mechanism  or  mechanisms  remain  unclear 
at  present.   It  is  the  intent  of  this  investigation,  within 
the  larger  context  of  the  Anomaly  Dynamics  Study  (ADS) ,  to 
"Describe  and  understand  the  mechanisms  responsible  for 
large  scale  climatological  anomalies  and  related  changes  in 
the  thermal  structure  and  the  general  circulation  of  the 
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upper  layers  of  the  ocean  in  order  to  predict  their  forma- 
tion and  evolution."  (Kirwan,  1976)   The  contribution  of 
anomalous  atmospheric  wind  forcing  as  a  generating  mechanism 
for  changes  in  the  thermal  structure  of  the  ocean  will  be 
emphasized . 

Previous  model  investigations  of  the  dynamics  of  large 
scale  temperature  anomaly  formation  in  the  ocean  have  been 
limited  by  the  lack  of  verifying  data,  at  the  surface  and 
at  depth,  over  sufficiently  large  areas.   In  cases  where 
anomalous  atmospheric  forcing  for  a  certain  period  of  time 
was  added  to  climatological  wind  forcing  (Namias,  1965; 
Arthur,  1966;  Jacob,  1967) ,  model  results  show  an  important 
influence  on  the  ocean  by  the  anomalous  forcing. 

In  the  present  study,  an  ocean  circulation  model  is 
used  to  make  a  number  of  initial  value  predictions  of  the 
evolution  of  large  scale  ocean  anomalies.   The  different 
model  predictions  are  designed  to  examine  the  effect  on  the 
prediction  of  including  or  excluding  certain  physical  pro- 
cesses in  the  model.   The  study  is  unique  in  that  a  care- 
fully collected  and  analyzed  set  of  ocean  temperature  data 
are  used. 

It  must  be  emphasized  at  this  point  that  this  study  is 
not  meant  to  be  a  test  of  the  predictive  capability  of  the 
model  itself.   First,  a  very  important  part  of  the  atmos- 
pheric forcing  fields,  namely  anomalous  heat  flux  that  is 
required  for  the  testing  of  the  prediction  capabilities  of 
any  ocean  model  is  not  available  on  the  time  or  space  scales 
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required.   Consequently,  only  anomalous  winds  are  used  in  the 
present  studies.   Second,  ocean  modeling  itself  is  still 
relatively  unsophisticated  and  certainly  untested  in  model- 
ing large  scale  thermal  anomalies  observed  in  the  oceans 
even  if  all  the  atmospheric  forcing  were  known.   In  addition, 
this  study  will  serve  as  an  initial  investigation  into  the 
adequacy  of  currently  available  temperature  data  in  describ- 
ing the  true  thermal  condition  of  the  Central  Pacific.   With 
the  current  oceanic  and  atmospheric  data  available  and  with 
the  current  state  of  the  science  of  ocean  modeling,  all  that 
can  be  tested  are  a  number  of  hypotheses  regarding  the  rela- 
tive importance  of  various  dynamic  forces  that  act  to  drive 
the  oceans.   Any  conclusions  resulting  from  this  study  then 
are  subject  to  re-evaluation  as  ocean  modeling  and  monitoring 
capabilities  improve.   The  1976-77  fall-winter  period  was 
selected  for  this  study  because  of  the  large  amplitude  anoma- 
lies in  the  atmosphere  and  ocean  that  existed  in  the  North 
Pacific  at  that  time. 
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II.    DATA  PROCESSING 

The  modeling  experiments  described  below  made  use  of 
anomalous  atmospheric  winds  and  oceanic  temperatures.   The 
first  step  in  preparing  the  atmospheric  data  for  input  into 
the  ocean  model  was  calculation  of  sea-level  pressure  anoma- 
lies or  departures  from  normal  for  all  of  1976  and  January 
and  February  of  1977.   All  the  atmospheric  data  were  generous- 
ly made  available  to  us  by  Dr.  J.  Namias  of  Scripps.   Long 
term  (20  year)  mean  values  of  sea  level  pressure  for  each  of 
these  months,  as  determined  by  Namias,  were  subtracted  from 
monthly  means  in  order  to  obtain  the  sea  level  pressure  de- 
partures from  norm.   Analysis  of  these  anomaly  fields  indi- 
cated the  period  from  September  1976  until  January  1977  ex- 
hibited sufficiently  large  anomalous  sea  level  pressures  to 
be  suitable  for  our  investigation. 

In  order  then  to  input  Namias'  anomalous  pressure  data 
into  the  model,  it  was  necessary  to  extrapolate  and  then 
interpolate  the  values  from  Namias'  5  °-  latitude  by  5°  longi- 
tude (9x19)  fields  encompassing  20°N  to  60°N  and  145°E  to 
125°W  to  the  33x33  fields  extending  from  the  equator  to  65°N 
and  145°E  to  125°W  required  by  the  model  (see  figure  1)  . 
The  extrapolation  of  the  anomalous  pressure  data  toward  the 
northern  boundary  assumed  a  slope  that  decreased  by  50  per- 
cent at  each  grid  point.   Extrapolation  towards  the  equator 
was  such  that  anomalous  u  and  v  components  of  the  wind  field, 
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calculated  from  the  anomalous  pressures,  go  to  zero  at  the 
equator.   Once  the  extrapolation  was  completed,  data  north 
of  50°N  was  smoothed  several  times  in  the  east-west  direc- 
tion in  order  to  remove  any  small  scale  features  introduced 
by  the  extrapolation.   This  procedure  was  carried  out  for 
the  months  of  September  1976  through  January  1977. 

With  the  33x33  field  of  sea  level  pressure  anomalies 
computed,  the  corresponding  anomalous  geostrophic  u  and  v 
components  of  the  wind  were  calculated.   These  anomalous  geo- 
strophic wind  fields  were  the  only  anomalous  atmospheric 
conditions  used  in  the  ocean  model.   These  fields  (figure  2) 
are  described  in  the  Results  section  below. 

The  next  step  in  data  preparation  for  input  into  the 
model  was  to  compute  the  temperature  anomalies  that  would 
serve  as  initial  and  verifying  conditions  for  the  model. 
All  temperature  data  was  provided  by  TRANSPAC,  the  ships  of 
opportunity  cooperative  research  effort.   The  TRANSPAC  ex- 
pendable bathythermograph  (XBT)  field  program  is  a  part  of 
the  North  Pacific  Experiment  (NORPAX)  which  is  charged  with 
investigating  how  oceanic  variability  influences  the  weather 
and  climate,  and  in  turn  how  the  weather  and  climate  varia- 
bility influence  the  oceanic  circulations.   Prior  to  TRANSPAC, 
data  on  the  thermal  structure  of  the  ocean  was  inadequate, 
thus  leading  to  the  inception  of  the  TRANSPAC  XBT  field  pro- 
gram to  provide  the  subsurface  temperature  data  necessary 
for  NORPAX.   The  observational  phase  of  the  program  began  in 
January  1975  utilizing  ships  of  opportunity  (i.e.  merchant 
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ships  voluntarily  providing  a  no-cost  observation  platform) 
that  routinely  traverse  the  sea  lanes  between  the  United 
States  and  Japan.   The  temperature  anomaly  data  provided 
through  TRANSPAC  was  obtained  by  subtracting  climatology 
means  from  the  monthly  mean  observed  temperatures.   Anomalous 
temperatures  were  provided  at  ten  levels  including  0,30,  60, 
90,  120,  150,  200,  250,  300  and  400  meters.   For  initializ- 
ing temperature  fields  in  the  model,  the  anomalous  tempera- 
ture field  for  September  1976  was  added  to  the  model  gener- 
ated climatology  for  September.   Calculation  of  verification 
statistics  entailed  use  of  anomalous  temperature  fields  for 
October  through  December  of  1976  and  January  1977.   Again, 
the  TRANSPAC  anomaly  data  had  to  be  processed  for  use  in  the 
model  because  the  domain  and  standard  levels  of  the  TRANSPAC 
data- did  not  coincide  with  the  domain  and  standard  levels  of 
the  model.   It  was  decided,  in  order  to  avoid  major  changes 
to  the  model,  to  extrapolate  from  the  2°  latitude  by  5°  longi- 
tude (11x15)  TRANSPAC  fields  extending  from  30°N  to  50°N 
and  160°E  to  130°W  to  the  area  and  gridpoints  encompassed 
by  the  model  (see  figure  1). 

The  first  step  in  preparing  the  TRANSPAC  temperature  anom- 
aly data  was  to  interpolate  the  data  from  the  ten  levels  of 
TRANSPAC  to  11  intermediate  levels  of  the  model,  specifically 
to  0,  20,  45,  80,  125,  200,  325,  600,  1200,  2200  and  4000 
meters.   Below  400  meters,  values  were  assigned  to  the  inter- 
mediate levels  by  decreasing  the  400  meter  anomalous  tempera- 
ture to  a  value  of  zero  at  the  2200  meter  intermediate  level. 
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This  was  accomplished  by  multiplying  the  400-  meter  value  by 
a  cosine  function  which  decreased  with  depth.  For  example, 
if  T1  (400)  is  the  given  TRANSPAC  anomaly  at  any  grid  point, 
then  the  anomaly  at  depth  Z,  T ' (Z) ,  was  specified  by 

T'(Z)     =    T'(400)     cos      [(Zi840°0°)     f]      • 

This  is  a  very  important  assumption,  since  it  assumes  the 
temperature    anomaly  cannot  change  sign  with  depth  below 
400  meters.   The  resulting  temperatures  at  the  11  intermedi- 
ate levels  were  then  averaged  between  any  two  levels  to  de- 
fine the  temperature  at  the  ten  prognostic  levels  in  the 
model  which  are  located  at  10.0,  32.5,  62.5,  102.5,  162.5, 
262.5,  462.5,  900.0,  1700.0  and  3100.0  meters. 

The  next  step  in  preparing  the  TRANSPAC  data  was  to  ex- 
pand it  to  cover  the  entire  area  of  the  model.   This  was 
done  by  expanding  the  temperature  anomaly  at  the  border  of 
the  11x15  TRANSPAC  area  outward  to  zero  at  the  outer  boundary 
of  the  model.   In  order  to  provide  a  smooth  transition  into 
the  area  of  no  data,  the  extrapolation  was  along  a  sine  func- 
tion as  shown  in  figure  1.   Areas  of  the  array  which  could 
be  filled  by  both  longitudinal  and  latitudinal  extrapolation 
were  assigned  a  temperature  anomaly  value  which  was  an 
average  of  the  two  extrapolations.   Once  the  entire  geographi- 
cal area  of  the  model  was  filled  on  a  2°  latitude  by  5° 
longitude  grid,  an  interpolation  scheme  was  used  to  fill 
the  33x33  array  of  model  grid  points.   The  33x33  array  of 
anomalous  temperatures  at  the  ten  levels  of  the  model  for 
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September  1976  was  then  added  to  the  model  generated  clima- 
tology for  September  15  to  constitute  the  initial  tempera- 
ture conditions  of  the  model.   Similar  33x33x10  arrays  for 
October  1976  to  January  1977  served  as  the  verifying  data 
for  computing  verification  statistics  for  the  different 
model  integrations  as  described  below. 


16 


Ill .    MODELING  APPROACH 

In  order  to  test  a  variety  of  atmospheric  and  oceanic 
processes  that  may  be  responsible  for  the  formation  and 
evolution  of  the  large-scale  temperature  anomalies  that  are 
of  interest,  a  series  of  initial  value  model  integrations 
have  been  carried  out  using  a  10-level  primitive  equation 
model  in  a  closed  rectangular  basin.   The  model  is  that  of 
Haney  (1974)  which  has  been  improved  to  include  time  de- 
pendent seasonal  forcing  by  the  atmosphere,  a  parameteriza- 
tion of  surface  wind  and  convective  mixing  (Haney  and 
Davies,  1976) ,  and  nonlinear  lateral  eddy  viscosity  based 
on  two-dimensional  turbulence  theory  (Haney  and  Wright, 
1975)  .   The  ocean  model  is  driven  by  wind  and  differential 
heating  which  is  calculated  from  the  model  sea-surface 
temperature  and  prescribed  values  of  the  atmospheric  solar 
radiation,  total  cloudiness,  surface  air  temperature,  rela- 
tive humidity  and  winds.   These  quantities  are  made  up  of 
a  climatological  part  and  in  some  cases  an  anomalous  part. 

To  initialize  the  model  for  a  prediction  experiment, 
the  model  is  first  integrated  over  a  long  "spin-up"  period 
(240  years)  using  climatological  atmospheric  forcing.   The 
last  year  of  this  integration  is  referred  to  as  the  "model 
climatology"  and  the  resulting  currents  and  temperatures, 
especially  in  the  mid-ocean  areas,  bear  a  close  resemblance 
to  those  in  the  North  Pacific  Ocean.   The  initial  conditions 
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for  the  prediction  experiments  are  then  obtained  by  combining 
a  given  anomaly,  whether  oceanic  or  atmospheric,  with  the 
model  climatology  for  the  same  parameter  for  the  appropriate 
time  of  the  year. 

Once  the  anomalous  atmospheric  forcing  and  temperature 
anomaly  fields  were  prepared  for  input  into  the  model,  five 
runs,  each  with  different  model  configurations,  were  made. 
As  stated  in  the  introduction,  the  purpose  of  the  different 
model  experiments  was  to  examine  the  relative  importance  of 
various  atmospheric  and  oceanic  processes  in  the  formation 
and  evolution  of  the  observed  anomalies. 

The  five  different  model  integrations  are  most  easily  de- 
scribed by  writing  the  equation  of  horizontal  motion  and  the 

thermodynamic  energy  equation  as  follows,  (symbols  defined  in 
Table  I) , 

dVH     1  3     8VH 

XT1  =  T  7»P  *  fK  x  vn  +  V«(AVV„)  +  -^(K  T— )        (1) 
dt      DH  H  H      dZdZ 

O 

dT  3     St 

~   =  M(Q,T)  +  V-  (AVT)   +  -^(K  -||)  (2) 

The  terms  on  the  right  hand  side  of  (1)  are  the  pressure 
gradient  force,  the  Coriolis  force,  horizontal  momentum  dif- 
fusion and  vertical  diffusion.   The  terms  on  the  right  hand 
side  of  (2)  are  the  vertical  mixing  of  heat  due  to  the  action 
of  surface  wind  and  convection,  horizontal  diffusion  and 
vertical  diffusion.   The  ocean  model  is  driven  by  a  downward 
flux  of  heat  and  momentum  (stress)  at  the  surface  which 
enters  through  boundary  conditions  in  the  vertical  diffusion 
terms : 
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(K  21) 

3Z  Z=0   =  Q  =  downward  heat  flux  at  the  surface 

3v 

H 
(K  , _  ) _    =  T  =  surface  stress 
a  Z   Z  =  0 


In  all  of  the  model  experiments  climatological  values  of  the 
heat  flux,  0  ,  were  used.   The  stress/  T  ,  was  calculated  from 


T  =  p  c:   V  I  V   , 
Ka  D '  S '   S 


where   V    is  a  prescribed  surface  wind  which  is  the  sum  of 
a  climatological  part  and  an  anomalous  part;  the  latter 
being  calculated  from  Namias '  sea  level  pressure  anomalies 
described  above.   Thus,  only   V    has  an  anomalous  part,  say 
V  '   and  it  enters  the  model  equations  in  two  places;  as  a 

boundary  condition  (wind  stress)  in  the  vertical  diffusion 

8VH 
term,  (K  t- — )      ;  and  directly  in  the  vertical  mixing  term 
dZ   Z  —  0 

M(Q,T) .   The  ability  exists  in  the  model  then  to  allow  the 

anomalous  surface  winds  to  enter  into  either,  both  or 

neither  of  the  terms.   In  this  way,  we  examine  separately 

the  role  of  anomalous  wind  driven  currents  and . anomalous 

wind  generated  vertical  mixing  in  the  generation  of  oceanic 

thermal  anomalies. 

In  physical  terms,  the  first  model  configuration,  model 

1,  describes  the  free  evolution  of  an  initial  temperature 

anomaly  since   V  '  =  0   in  this  model.   In  model  2,  the 

S 

effect  of  the  anomalous  wind  stress  in  driving  anomalous 
Ekman  type  flow  is  added,  since   V  '  =  0   only  in  the  mixing 
term  in  this  model.   According  to  steady  Ekman  theory  this 
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flow  would  be  90°  to  the  right  of  the  anomalous  surface 
winds.   Model  3  incorporates  the  effects  of  anomalous  wind 
stress  in  both  Ekman  type  flow  and  surface  wind  mixing. 
This  mixing  process  is  expected  to  produce  increased  verti- 
cal mixing  in  regions  of  strong  anomalous  winds.   Model  4 
is  similar  in  its  dynamics  to  model  2  but  has  the  initial 
anomalous  temperature  fields  equal  to  zero.   Likewise, 
model  5  is  similar  to  model  3  but  with  the  initial  anomalous 
temperatures  equal  to  zero.   These  last  two  experiments  are 
designed  to  evaluate  the  importance  of  having  good  initial 
conditions  when  making  a  prediction. 

In  order  to  provide  a  quantitative  comparison  between 
the  different  models,  a  variety  of  statistics  were  computed 
for  the  five  model  cases  at  30-day  intervals  of  time.   The 
statistics  calculated  included  the  average  absolute  predic- 
tion error,  root  mean  square  error  and  an  SI  score  associated 
with  the  model  generated  temperature  anomalies.   The  follow- 
ing definitions  indicate  exactly  what  was  measured  by  the 
statistics.   Average  absolute  error  and  root  mean  square 
error  are  defined  as  follows:   if  in  a  series  of  N  forecasts 
(i.e.  N  different  gridpoints),  F.   represents  the  i-th  fore- 
cast and   0.   the  corresponding  observation,  the  average 
absolute  error  is  given  by 


Z  F.  -  0. 

1  i i_ 

N 
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and  the  root  mean  square  error  by 


RMSE  = 


V 


2(Fi  -  o±) 

N 


(Brier  and  Allen,  1951)  . 

The  SI  score  is  a  dimens ionless  measure  of  the  difference 
between  predicted  and  observed  horizontal  gradients  and  is 
given  by 

EleJ 


SI  = 


Z    e 


(Teweles  and  Wobus,  1954), 

where   e    is  the  error  in  the  forecast  temperature  anomaly 
G  tr  j 

gradient  in  the  horizontal  and   G    is  the  observed  tempera- 

L 

ture  anomaly  gradient.   The  SI  score  was  further  divided  in 
order  to  measure  the  skill  in  forecasting  both  the  north- 
south  and  the  east-west  gradient  of  the  temperature  anomaly. 

The  various  model  predicted  and  TRANSPAC  derived  tem- 
perature anomaly  fields  were  compared  beginning  with  the 
October  1976  fields  at  30-day  intervals  of  time  and  at  nine 
levels  of  the  model.   Statistics  were  computed  only  over  an 
area  from  35°N  to  45°N  and  130°W  to  180°.   The  verification 
was  restricted  to  the  region  east  of  the  dateline  because 
the  region  west  of  the  dateline  is  generally  characterized 
by  numerous  mesoscale  temperature  features  (Bernstein  and 
White,  1977).   Since  the  intent  of  this  investigation  was  in 
part  to  study  present  capabilities  in  modeling  large  scale 
thermal  anomalies,  it  was  assumed  the  model  could  not  effec- 
tively predict  such  mesoscale  features. 
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IV.    RESULTS  AND  CONCLUSIONS 

A  vast  amount  of  data  is  made  available  by  this  study 
and  as  is  always  the  case,  interpreting  such  a  volume  of 
data  in  light  of  the  original  intentions  of  the  study  is 
difficult.   Describing  the  actual  evolution  of  the  thermal 
anomalies  in  the  observed  and  various  model  cases  is  perhaps 
the  most  logical  starting  point.   A  later  more  detailed 
treatment  of  the  verification  statistics  themselves  will 
serve  to  highlight  various  aspects  of  the  described  results. 

Since  the  effects  of  anomalous  wind  stress  in  producing 
large  scale  thermal  anomalies  is  considered  paramount,  a 
description  of  the  anomalous  atmospheric  forcing  for 
September  1976  to  January  1977  should  precede  any  descrip- 
tion of  observed  or  model  evolutions  of  thermal  anomalies. 
Figure  2  shows  the  anomalous  sea-level  pressure  fields  for 
these  months,  contoured  at  2  mb  intervals  with  shaded  areas 
indicating  negative  anomalous  pressures.   These  anomalous 
pressure  fields  are  shown  only  over  the  area  of  the  tem- 
perature fields  given  in  figures  3  through  8.   The  September 
1976  mean  circulation  over  the  central  Pacific  was  dominated 
by  a  trough  at  700  mb  which  had  progressed  to  a  position 
south  of  the  Aleutians.   Downstream,  the  western  North 
American  ridge  became  established  as  an  Arctic  ridge 
(Taubensee,  1976).   This  mean  700  mb  trough  in  the  central 
Pacific  is  evidenced  at  the  surface  by  negative  anomalous 
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pressure.   In  October  1976,  the  mean  700  mb  trough  over  the 
eastern  Pacific  retrograded  somewhat,  as  did  the  downstream 
ridge  over  western  North  America  (Wagner,  1977a).  The  nega- 
tive anomalous  pressure  pattern  at  sea  level    is  seen  in 
figure  2  to  move  westward  in  October,  being  replaced  to  the 
east  by  positive  anomalous  pressure  associated  with  the 
ridge  over  western  North  America.   The  mean  700  mb  trough 
over  the  central  Pacific  remained  relatively  stationary  but 
intensified  significantly  (2  standard  deviations  from  the 
mean)  from  November  1976  to  January  1977  (Dickson,  1977; 
Taubensee,  1977;  Wagner,  1977b)  .  This  intensification  is  mani 
fested  in  the  sea  level    pressure  by  an  extreme  deepening 
of  the  negative  anomalous  pressure  fields  for  November 
through  January. 

The  observed  anomalous  temperature  data  for  September 
1976,  November  1976  and  January  1977  at  10,  102,  162  and  262 
meters  are  shown  in  figure  3.   The  observed  evolution  during 
this  winter  period  in  the  upper  level  is  characterized  by 
extreme  cooling  over  all  but  the  easternmost  edge  of  the 
area.   In  the  mid-levels  (102  and  162  meters) ,  the  September 
field  shows  an  extremely  intense  disturbance  in  the  monthly 
mean  temperature  anomaly  at  and  just  east  of  the  dateline. 
This  disturbance  propagates  westward  out  of  the  area  by 
November.   The  102  meter  level  exhibits  a  general  cooling 
over  the  winter  months,  similar  to  that  in  the  upper  level. 
At  162  meters  there  is  a  general  trend  towards  warming  with 
the  axis  of  temperature  anomalies  shifting  from  generally 
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north-south  to  east-west.   The  262  meter  field  shows  a  less 
intense  anomaly  structure  than  the  mid-levels  for  September 
but  exhibits  a  more  intense  warming  along  the  east-west  axis 
than  had  occurred  by  January  at  162  meters. 

Figure  4  shows  the  anomaly  evolution  predicted  by  model  1 
(climatological  winds) .   The  dominant  feature  in  the  upper 
levels  is  the  development  of  an  intense,  basically  mesoscale 
feature  in  the  western  portion  of  the  area.   This  development 
is  due  to  the  strong  north-south  advection  of  mean  tempera- 
ture by  the  anomalous  geostrophic  surface  currents  induced 
in  the  model  by  the  intense  disturbance  in  the  western  part 
of  the  basin  at  mid-levels.   This  same  type  of  development 
occurs  in  the  middle  and  lower  levels  so  that  a  similar  but 
less  intense  structure  develops  in  the  vertical.   These 
features  tend  toward  a  slow  westward  propagation  during  the 
winter  months,  indicating  partial  development  and  movement 
as  free  Rossby  waves. 

Model  2  introduces  anomalous  winds,  but  only  in  the  sur- 
face stress  term  (climatological  winds  are  used  in  the  verti- 
cal mixing  term) .   Including  this  effect  produces  anomalous 
Ekman  type  flow  in  the  upper  levels  which  definitely  improves 
the  model  predictions  by  cooling  the  upper  layers  (see  figure 
5)  .   The  intense  mesoscale  features  in  the  west  produced  by 
the  strong  north-south  advection  of  mean  temperature  de- 
scribed in  model  1  are  still  very  apparent  in  the  anomaly 
structure  however.   These  anomalously  cold  mesoscale  features 
are  intensified  and  expanded  in  their  horizontal  scale  in  the 
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upper  levels  by  the  superposition  of  the  large  scale  cooling 
resulting  from  the  anomalous  wind  stress.   In  the  middle  to 
lower  levels  there  is  a  general  warming,  similar  to  the  ob- 
served case,  which  tends  to  reduce  the  intensity  of  the  cold 
eddies . 

Model  3  includes  the  effects  of  anomalous  winds  in  both 
the  surface  stress  term  and  the  vertical  wind  mixing  of  heat. 
This  wind  mixing  process  would  be  expected  to  produce  in- 
creased vertical  mixing  in  regions  of  increased  surface  winds 
but  as  seen  in  figure  6,  the  results  for  model  3  do  not  vary 
significantly  from  model  2.   In  the  ocean  model,  surface  wind 
mixing  depends  in  part  on  the  net  surface  heat  flux.   If  the 
net  flux  is  downward,  wind  mixing  occurs,  but  if  the  net  sur- 
face heat  flux  is  upward  as  it  is  in  this  winter  case  there 
is  little  or  no  wind  mixing.   Model  3  results  then  show  only 
slight  differences  from  model  2  in  the  uppermost  layers. 

Model  4  includes  anomalous  winds  only  in  the  surface 
stress  term  (like  model  2)  but  the  initial  temperature  anomaly 
is  set  to  zero  in  order  to  evaluate  the  importance  of  knowing 
the  initial  conditions.   In  this  particular  winter  case  of 
very  strong  anomalous  wind  forcing,  the  temperature  anomalies 
in  the  upper  levels  are  predicted  slightly  better  than  in 
models  2  or  3,  with  general  cooling  over  most  of  the  area  as 
seen  in  figure  7.   Results  in  middle  to  lower  levels  are 
worse,  however,  with  the  102  m  level  showing  development  of 
small  scale  warm  anomalies  near  30-35N  where  cooling  actually 
occurred  and  the  262  m  level  showing  cooling  to  the  northeast 
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where  warming  is  observed.   It  may  be  noted  at  this  point 
that  the  superposition  of  the  results  for  model  4  with 
initial  anomalous  temperature  equal  zero  and  model  1  with 
no  anomalous  wind  forcing  will  give  almost  exactly  the 
results  observed  in  model  2.   This  is  cle"ar  evidence  that 
in  the  model,  the  effects  of  external  (wind)  forcing  and  the 
effects  of  internal  dynamical  development  take  place  rather 
independent  of  one  another. 

Model  5  is  the  same  in  its  dynamics  as  model  3  but  again 
the  initial  anomalous  temperature  is  set  to  zero.   As  pre- 
viously discussed,  models  2  and  3  show  little  variation  and 
as  seen  in  figure  8,  model  5  results  vary  only  slightly  from 
model  4. 

The  actual  verification  statistics, of  course, provide  more 
concrete  evidence  to  the  similarity  and  differences  between 
the  various  models  in  successfully  describing  the  evolution 
of  the  large  scale  temperature  anomalies  as  observed.   Table 
TI  presents  these  statistics  in  the  following  format:  nine 
levels  from  10  meters  to  1700  meters  depth  are  considered; 
the  ratio  of  the  root  mean  square  error  (RMSE)  to  the  RMSE 
of ' climatology  (a  prediction  of  zero  temperature  anomaly 
everywhere)  is  given  as  the  first  number  followed  by  a 
similarly  non-dimensionalized  SI  score  for  30  day  intervals 
from  October  1976  to  January  1977.   A  ratio  of  RMSE ' s  and  Si's 
were  used  as  a  measure  of  skill  for  easy  interpretation  since 
a  ratio  less  than  one  indicates  skill  over  climatology 
whereas  a  ratio  greater  than  one  indicates  less  skill  than 
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climatology.   The  skill  of  a  forecast  of  persistence  is  also 
shown.   These  ratios  for  10,  102,  162,  and  262  meters  have 
been  plotted  as  a  function  of  time  in  figures  9a  and  9b  for 
easy  reference. 

The  plot  of  skill  as  a  function  of  time  at  10  meters 
seen  in  figure  9a  indicates  an  extremely  rapid  loss  of  skill 
as  early  as  30  days,  followed  by  improved  skill  in  the  models 
(2,  3,  4  and  5)  having  wind  forcing.   Model  1  with  no 
anomalous  wind  forcing  shows  no  skill  after  30  days  but 
models  2  and  3  clearly  demonstrate  the  benefits  of  consider- 
ing such  anomalous  forcing.   Surprisingly,  models  4  and  5 
with  zero  initial  anomalous  temperature  conditions  show  con- 
stant improvement  with  time  even  over  models  2  and  3.   This 
would  apparently  indicate  the  greater  importance  of  knowing 
the  anomalous  wind  forcing  than  of  knowing  the  initial  tem- 
perature field  exactly,  especially  in  a  case  as  this  one 
where  the  anomalous  winds  were  very  strong. 

Perhaps  the  most  conspicuous  result  is  the  fact  that  the 
model  predictions  at  the  surface  are  better  when  the  initial 
anomaly  is  zero  (model  4  and  5)  than  when  it  is  the  observed 
anomaly  (model  2  and  3) .   In  the  latter  case,  strong  geo- 
strophic  flow  in  the  north-south  direction  is  induced  at 
the  surface  and  this  flow  advects  the  mean  isotherms  in  the 
model  (cold  to  the  north  and  warm  to  the  south)  to  produce 
a  large  amplitude  disturbance  in  the  western  third  of  the 
basin  which  does  not  exist  in  the  observed  pattern.   As  de- 
scribed earlier,  the  strong  meridional  flow  is  due  to  the 
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intense  anomaly  in  the  initial  temperature  field  at  and  be- 
low 102  meters.   One  possible  explanation  for  the  model 
discrepancy  is  that  the  model  simply  does  not  propagate  this 
developing  disturbance  westward  and  out  of  the  domain  fast 
enough.   However,  an  examination  of  the  surface  anomaly 
data  west  of  the  dateline  shows  no  large  disturbance  there 
to  which  the  predicted  disturbance  in  figure  4  could  corre- 
spond if  the  model  error  was  simply  one  of  phase  error. 
A  second  possibility  is  that  in  September  there  existed  a 
compensating  deep  temperature  anomaly  (below  400  meters) 
whose  induced  vertical  shear  nearly  cancels  that  induced  by 
the  anomaly  observed  between  102  and  262  meters.   In  this 
case,  the  model  assumption  that  the  initial  anomaly  decreases 
slowly  to  zero  at  depth,  without  changing  sign  (see  DATA 
PROCESSING)  is  wrong.   If  this  second  explanation  is  true, 
it  implies  that  in  order  to  predict  or  understand  sea  surface 
temperature  anomalies,  one  must  measure  the  anomaly  field 
all  the  way  down  to  the  bottom  of  the  sea.   Unfortunately, 
it  will  require  many  case  studies  of  this  type  before  this 
question  is  adequately  resolved. 

The  plot  of  skill  versus  time  for  102  meters  in  figure 
9a  shows  a  somewhat  less  rapid  loss  of  skill  initially  but 
demonstrates  a  fairly  high  amplitude  oscillation  about  1.0 
in  the  30  to  120  day  period  of  time.  Comparison  of  the 
curves  for  models  4  and  5  and  models  1,  2  and  3  indicates 
a  relation  between  these  oscillations  and  the  initial  tem- 
perature conditions  since  this  alternating  behavior  is  not 
displayed  by  models  4  and  5. 
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Figure  9b  shows  an  extreme  and  rapid  loss  of  skill  (com- 
pared to  climatology  i.e.  a  forecast  of  zero  anomaly)  for 
persistence  and  the  models  at  162  and  262  meters.   An  oscil- 
lation in  performance  characterizes  the  162  meter  graph. 
It  is  important  to  note  also  that  while  models  4  and  5  have 
greater  skill,  they  are  not  as  good  as  climatology.   The  262 
meter  plot  does  not  exhibit  this  oscillatory  characteristic, 
presumably  as  a  result  of  less  intense  anomalous  initial 
temperatures  at  this  depth.   It  is  rather  disheartening  that 
the  model  loses  skill  so  fast  at  these  levels.   It  is  sur- 
prising too  that  persistence  is  also  a  very  bad  forecast 
at  these  depths  which  are  generally  believed  to  undergo 
slow  changes. 

The  SI  scores  computed  to  measure  skill  in  predicting  the 
gradient  of  large  scale  thermal  anomalies  are  shown  in  Table 
I.   Models  4  and  5  consistently  demonstrate  greater  skill 
than  either  persistence  or  the  other  models  in  predicting 
gradients.   This  indicates  of  course  that  the  ocean  model  has 
little  skill  in  computing  gradients  since  any  field  more 
closely  resembling  a  zero  field  possesses  greater  skill.   If 
the  model  SI  scores  are  broken  up  into  north-south  and  east- 
west  components,  however,  there  is  greater  relative  skill  in 
predicting  north-south  gradients  than  east-west  gradients. 

The  results  of  this  study  can  be  summarized  as  follows. 
In  the  upper  levels,  the  models  which  include  anomalous  wind 
forcing  show  skill  over  climatology  out  to  4  months,  i.e.  the 
ratio  of  RMSE  of  the  model  to  RMSE  of  climatology  is  less 
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than  one.   However,  knowledge  of  the  initial  anomalous  tem- 
perature conditions  degrades  the  prediction  skill.   Two 
possible  explanations  were  offered  but  a  variety  of  addi- 
tional case  studies  are  needed  to  properly  resolve  this 
problem.   In  the  lower  layers  (162  to  262  meters),  none  of 
the  models  show  skill  over  climatology  and  consideration 
of  neither  anomalous  wind  forcing  nor  initial  anomalous  tem- 
peratures improves  the  models  predictive  capabilities.   The 
rather  poor  model  results  at  these  depths  are  a  serious  con- 
cern because  it  is  not  likely  that  the  addition  of  anomalous 
heat  fluxes  to  the  model  predictions  will  have  much  effect 
at  these  levels.   In  the  future,  it  is  planned  to  repeat 
this  type  of  study  for  cases  in  which  the  anomalous  atmo- 
spheric forcing  was  very  weak  and  also  to  repeat  the  present 
study  adding  anomalous  heat  fluxes  being  computed  by  other 
NORPAX  investigations. 
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V  =  horizontal  velocity 
H 

p  =  reference  density 
o 

p  =  pressure 

f  =  Coriolis  parameter 

A  =  horizontal  diffusion  coefficient 

K  =  vertical  diffusion  coefficient 

T  =  temperature 

M  =  parameterized  vertical  mixing  term 

Q  =  surface  heat  flux 

T  =  surface  wind  stress 

TABLE  I 
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10  M  (RMSE   RATIO/SI  RATIO) 

OCT               N0V               DEC             JAN 

PERS 

1.17/1.53 

1. 11/1.42 

1.12/1.61 

1.18/1.76 

Ml 

1.12/1.46 

1.08/1.61 

1.07/1.77 

1.11/1.79 

M2 

1.04/1.48 

0.89/1. 58 

0.70/1. 56 

0.77/1. 52 

M3 

1.06/1.53 

0.87/1.66 

0.68/1.57 

0.77/1. 51 

M4 

0.92/0.99 

0.82/1.05 

0.63/0.99 

0.69/1.01 

M5 

0.96/1.13 

0.80/1.14 

0.62/0.99 

0.69/0.99 

32  M  (RMSE  RATIO/SI  RATIO) 

OCT              NOV              DEC             JAN 

PERS 

1. 18/1.53 

1.12/1.50 

1.12/1.47 

1.13/1.64 

HI 

1.14/1.48 

1 .06/1.68 

1.07/1.77 

1.08/1.85 

M2 

1.16/1.63 

1.01/1.82 

0.81/1.56 

0.83/1.55 

M3 

1.08/1.91 

0.98/1.90 

0.79/1.57 

0.83/1. 54 

M4 

1.02/1.08 

0.96/1. 38 

0.75/0.99 

0.77/0.95 

115 

0.93/1.43 

0.93/1.45 

0.73/0.98 

0. 77/0.93 

62  M  (RMSE  RATIO/SI  RATIO) 

OCT               NOV               DEC             JAN 

PERS 

1.35/1. 23 

0.97/1.80 

0.94/1.43 

0.99/1.83 

141 

0.96/1.15 

0.99/1.70 

0.98/1.56 

1.08/1.91 

M2 

0.99/1.21 

1.05/1. 94 

0.94/1. 75 

0.81/1.79 

M3 

1.04/1. 36 

1.07/2.01 

0.95/1.  74 

0.81/1.79 

M4 

1.01/0. 97 

1.05/1. 29 

0.98/1. 36 

0.74/1.20 

115 

1.06/1.09 

1.07/1.41 

0.99/1.  36 

0.74/1. 21 

Table    II 
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102  M  (RMSE  RATIO/SI  RATIO) 

OCT              NOV              DEC             JAW 

PERS 

0. 87/1. 22 

1. 38/1.67 

0.84/1.63 

0.95/2.29 

Ml 

0.91/1.17 

1.14/1. 52 

0.79/1.32 

1 .03/1.94 

M2 

0.91/1.18 

1.17/1.55 

0.86/1.45 

0.98/2.01 

M3 

0.94/1. 24 

1.24/1.58 

0.91/1.50 

0.98/2 .02 

M4 

1.00/0.99 

1. 03/1.05 

1.05/1. 21 

0.97/1.31 

M5 

1.06/1.10 

|l. 17/1. 14 

1.09/1. 26 

|  0.98/1.32 

162  M  (RMSE  RATIO/SI  RATIO) 

OCT               NOV               DEC             JAN 

PERS 

1.04/1.46 

1.85/2.07 

1.26/1. 80 

1.88/1.77 

Ml 

1.00/1. 31 

1.54/1.91 

0.91/1.34 

1. 59/1*66 

M2 

1.00/1.32 

1.46/1.89 

0.91/1.53 

1.65/1.74 

M3 

1.00/1. 32 

1.46/1.89 

0.91/1.53 

1.65/1.74 

M4 

1.00/1.03 

1.00/1. 10 

1.00/1. 15 

1.06/1.18 

M5 

1.00/1.03 

1.00/1.10 

J  1.00/1. 15 

1.06/1.18 

262  M  (RMSE RATIO/SI  RATIO) 

OCT               NOV               DEC             JAN 

PERS 

1.64/1.76 

1.73/2.29 

1.69/2.10 

0.97/1.38 

Ml 

1.57/2.03 

1. 55/2.11 

1.38/1.85 

1.03/1.38 

M2 

1.50/2.02 

1.45/2.09 

1. 31/2 .00 

1.05/1.45 

M3 

1. 50/2.02 

L. 45/2 .09 

1. 31/2.00 

1.05/1. 45 

M4 

1.00/1.02 

L. 00/1. 13 

1.00/1.18 

1.03/1. 11 

M5 

1.00/1.02 

1.00/1.13 

1.00/1.18 

1.03/1.11 

Table    Ii(cont) 
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462  M  (RMSE RATIO/SI  RATIO) 

OCT               NOV               DEC             JAN 

PERS 

1.62/1.87 

1.00/1.66 

1.06/1.81 

0.89/1.37 

Ml 

1.62/1.78 

0.88/1. 35 

0.94/1.38 

0.86/1.17 

M2 

1. 54/1. 77 

0. 88/1.  36 

0.89/1.41 

0.89/1.  24 

M3 

1.54/1.77 

0.88/1.36 

0.89/1.41 

0.89/1.24 

M4 

1.00/0 .98 

1.00/1.03 

1.00/1.07 

1.06/1.11 

M5 

1. 00/0. 98 

1.00/1.03 
1 

1.00/1.07 

1.06/1.11 

900  M  (RMSE  RATIO/SI  RATIO) 

OCT              NOV               DEC             JAN 

PERS 

1.45/1.68 

0. 93/1.44 

1.00/1.64 

0.86/1.31 

Ml 

1.45/1.58 

0. 79/1.19 

0.93/1.29 

0.86/1. 11 

M2 

1.45/1.56 

0.79/1.18 

0.93/1.30 

0 .86/1. 16 

M3 

1.45/1.56 

0.79/1.18 

0.93/1.30 

0.86/1.16 

M4 

1.00/0.99 

1. 00/1. 02 

1.00/1.06 

1.04/1.08 

M5 

1.00/0. 99 

1.00/1.02 

1.0  7/1.06 

1. 04/1.08 

1700  M  (RMSE  RATIO/SI  RATIO) 

OCT               NOV               DEC             JAN 

PERS 

1. 50/1.68 

0.80/1.44 

1.00/1.64 

0.82/1.31 

Ml 

1. 50/1.58 

0.80/1.18 

0.83/1.29 

0.82/1. 11 

M2 

1. 50/1.57 

0.80/1.17 

0.83/1.30 

0.82/1. 13 

M3 

1. 50/1. 57 

0.80/1.17 

0.83/1.  30 

0.82/1.13 

M4 

1.00/0.99 

1.00/0.99 

1.00/0.99 

1.00/1. 02 

M5 

1. 00/0.99 

1.00/0.99 

1 .00/0.99 

1. 00/1. 02 

Table    II (cont) 
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